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We study the vertical structure of the transition layer between an accretion disc 
and corona in the context of the existence of two-phase medium in thermally unstable 
regions. The disc is illuminated by hard X-ray radiation and satisfies the condition 
of hydrostatic equilibrium. We take into account the energy exchange between hot 
corona (~ 10 8 K) and cool disc (~ 10 4 K) through the radiative processes and due 
to thermal conduction. We make local stability analysis of the case with conductivity 
and we conclude that thermal conduction does not suppress thermal instability. In 
spite of continuous temperature profile T(t) there are regions with strong tempera- 
ture gradient where spontaneous perturbations can lead to cloud condensation in the 
transition layer. We determine the minimum size Xtc of such a perturbation. 
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1 INTRODUCTION 

The X-ray emission from active galactic nuclei (AGN) has 
been observed for more than thirty years. Although a lot of 
research were done on this topic in the case of radio-quiet 
AGN we still don't know neither the mechanism leading to 
formation of the hot X-ray emitting plasma nor the geome- 
try of the X-ray source (for a review, see Mushotzky, Done 
& Pounds 1993, Czerny 1994). 

Recent X-ray data show that a typical X-ray spectrum 
of a Seyfert 1 galaxy consists of a power law component with 
energy index ~ 0.8—1.0, high energy cutoff above ~ 200 keV, 
and a reflection component characterized by low-energy cut- 
off due to the absorption by cold matter, high-energy cut-off 
due to the Klein-Nishina effect and the presence of the Kq 
line originating in a cold medium (e.g. Pounds et al.1990, 
Zdziarski et al. 1995). 

Underlying power-low corresponds to emission from hot 
region and this radiation is most probably thermal (an- 
nihilation line for par production wasn't observed in any 
spectrum, Maciolek-Niedzwiecki et al. 1995). Second com- 
ponent describes reflection of the radiation of hot plasma 
by cold matter. Approximately half of the primary emis- 
sion is observed unprocessed and half of it is reflected. Fast 
variability of the line component observed in a number of 
sources (Yaqoob et al. 1996) shows that the cold medium 
is located cospatially with the hot one although some con- 



tribution from the reflection by a distant dusty /molecular 
torus (Krolik et al. 1994, Iwasawa et al. 1997) also might be 
present. 

An attractive scenario of the coexistence of the hot and 
cold medium is an optically thin corona above an optically 
thick accretion disc. The existence of such a corona was sug- 
gested by Liang & Price (1977). 

Main question about coronae refers to the manner how 
the energy dissipation proceeds in order to heat up the 
plasma to an (electron) temperature ~ 10 9 K. Usually the 
fraction of energy dissipated in coronae / is a free parame- 
ter of a model (Haard & Maraschi 1991 1993, Kusunose & 
Mineshige 1994, Svensson & Zdziarski 1994). Another ap- 
proach is to assume that the flux in the corona is generated 
by viscosity, similarly like in accretion disc. Parameteriza- 
tion by viscosity coefficient a allows to predict the ratio of 
the energy generated in the disc to that generated in corona. 
Nevertheless in this models full disc/corona vertical struc- 
ture is averaged by two hot and cold layers, and the place 
of transition is either adopted arbitrarily or based on crude 
approximations (Nakamura & Osaki 1993, Zycki et al. 1995, 
Witt et al. 1997). 

It is shown (Krolik, Mckee & Tarter 1981 hereafter 
KMT) that the illumination of the cold matter (~ 10 4 K) 
by hard X-rays leads to spontaneous stratification into hot 
and cold part due to the thermal instability. The two layers 
are in pressure equilibrium and the discontinuity in the tem- 
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perature and density reflects the discontinuity in the cooling 
mechanism: hot layer is cooled by Compton scattering while 
cold layer is cooled by atomic processes. This phenomenon 
underlies the spontaneous division of the flow into disc and 
a corona and therefore it is important to study it in some 
detail. 

The nature of this phenomenon is well understood due 
to simple analytical studies (KMT 1981, Begelman, McKee 
& Shields 1983) based on local determination of the temper- 
ature. Later studies which included the radiative transfer in 
the irradiated slab were not suitable for the purpose to an- 
alyze the problem since either the instability showed up in 
these codes as numerical problems (Raymond 1993) or the 
iteration had to be started below the instability zone (Sincell 
i Krolik 1997). Therefore we follow the basic semi - analyt- 
ical approach to the description of the effects of irradiation 
but we enrich the physics involved in the process. 

The vertical structure of such a transition between the 
Compton cooled and atomic cooled medium in the context 
of accretion disc/corona system was studied in our previous 
paper (Rozanska & Czerny 1996 hereafter RC96). The model 
was parameterized by the hard X-ray flux irradiating the 
transition layer from above and by the soft UV flux (thermal 
radiation of the cold disc) from below. We took into account 
the hydrostatic equilibrium and we assumed purely radiative 
energy exchange between the corona and the disc. 

The result of that paper suggested that as soon as there 
was enough of the hard X-ray emission liberated e.g. in the 
upper layer of the corona there was a well defined transition 
layer between the disc and the corona. In this zone two- 
phase medium may form, i.e cool clouds can coexist with 
hot plasma under the constant pressure. 

In the present paper we discuss the transition layer con- 
sidering also the effect of thermal conduction. For the first 
time this effect was studied in relation to dissipative corona 
by Maciolek-Niedzwiecki, Krolik & Zdziarski (1997). They 
studied the influence of the heat flux transport on the struc- 
ture of the corona and the transition between the hottest 
layers (10 9 K) and Compton heated disc atmosphere (10 8 K). 
Their approach did not allow to determine properly the ba- 
sis of the corona as they did not include atomic processes 
and were unable to reach cool disc layers cooled predom- 
inantly by lines and bound-free transitions. Our paper is 
therefore complementary to theirs as we describe more deep 
parts of the transition starting from Compton heated layers 
and reaching down the cool disc (10 4 K). 

We show that the thermal conduction allows to find very 
sharp but nevertheless continuous solution for density and 
temperature profile throughout a transition region, without 
the need for the two-phase medium. However, we study the 
local stability of such a solution and we show that this solu- 
tion is not thermally stable which supports the view of the 
existence of the two-phase medium in the transition zone, 
as suggested by RC96. 



2 ASSUMPTIONS OF THE VERTICAL 
STRUCTURE 

We consider the optically thick and geometrically thin ac- 
cretion disc illuminated on the top by hard radiation flux. 
The origin of this hard radiation in radio quiet objects is still 



unknown and usually two models are under debate. High en- 
ergy photons can rise through the Comptonization: (1) on 
thermal electrons in optically thin corona above the disc or 
in the innermost part of the optically thin disc (2) on non- 
thermal electrons in the corona or at the basis of the jet- like 
(albeit not highly collimated) flow. 

We do not specify which interpretation is the correct 
one. Instead, we treat the produced hard X-ray flux as ex- 
ternal at a given radius. Therefore, if the hard X-ray emis- 
sion is produced in the corona it means that we study the 
boundary layers between the upper, very hot coronal zone 
and cool disc. On the other hand, if the X-ray emission is 
generated somewhere in the innermost part of the accretion 
flow and illuminates the disc at larger radii we actually cal- 
culate the entire corona with the transition zone. To simplify 
the terminology we just call the studied layer 'the transition 
zone' in both cases. 

The incident flux of hard radiation F^ardo illuminating 
the upper surface of the transition layer is a free parameter 
in our calculations. It is only that part of whole radiation 
flux generated in X-ray source, which is directed towards 
the disc. We take into account the decrease of this flux by 
'absorption on the spot' which means full ekstinction of flux 
through real absorpion and through scattering: 

Fhard{r) = FhardO CXp( — t), (1) 

where r is the optical depth measured from the top of the 
slab downwards to the cold disc. We don't take into ac- 
count ekstinction of flux after one scatter, so the matter 
on considered optical depth does not become the source of 
radiation. Such an assumption is valid as long as absorp- 
tion is the dominant source of opacity. Recent papers on 
X-ray spectrum reflected from accretion disc suggest that 
in case of illumination almost 90% of incident luminosity 
is deposited within the cold disc, while the remaining 10% 
is reflected ( George & Fabian 1991, Matt, Perola & Piro 
1991, Haardt & Maraschi 1993). There are strong observa- 
tional evidences that absorption dominates in the illumi- 
nated matter (Magdziarz & Zdziarski 1995). So most of the 
X-ray energy is converted into soft radiation and reemited 
according to the formula: 

dF softl (r) = F hard (l - a)dr (2) 

We assume the albedo of hard radiation from the disc atmo- 
sphere to be constant with optical depth and equal a = 0.15. 

Another source of energy is the gravitational energy of 
the accreting matter itself liberated via the viscous dissipa- 
tion (Czerny & King 1989ab). The soft radiation flux gen- 
erated in this way increases upward throughout the disc ac- 
cording to the formula: 

dF som (r) = -IS^-PdT, (3) 

2 pKtot 

where fi is the angular velocity of keplerian motion, a is 
the viscosity parameter, as introduced by Shakura & Sun- 
yaev (1973), P means the gas pressure, and K tot is the total 
opacity coefficient. 

Total soft flux emitted by the disc is a sum of two com- 
ponents described above: 

dF soft (r) = dF softl (r) + dF soJt2 (r) (4) 
The spectrum of the soft radiation is that of a black 
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body with the temperature Tt,b whilst hard X-ray emission is 
assumed to be a power law with the energy index ole = —0.9 
(here F v ~ v aE ). We neglect the change of the spectral 
shape in any of the two components with the optical depth. 

Strong illumination by X-rays heats up the outer parts 
of the disc and the very hot, optically thin slab is created 
above the disc. Such a corona generally reduces the radiation 
pressure, so we assume for the transition zone the following 
equation of state: 



P = 



-pT, 



(5) 



with the value of molecular weight p, — 0.5 for cosmic 
chemical composition, i.e. we neglect the radiation pressure 
throughout the zone. 

The vertical distribution of the pressure and the density 
p is defined by the condition of the hydrostatic equilibrium 
in vertical direction: 



p dz 



- "oft , 



-K-tot 



(6) 



Ha is the half of the disc thickness and the coordinate z is 
measured from the disc half-thickness, Hd, downwards. 

The temperature of a slab is determined by balanc- 
ing radiative heating and cooling and heat transport via 
the thermal conduction from the neighbor layers. Time- 
independent energy balance equation in vertical direction 
becomes: 



(7) 



where C is generalized loss-heat function defined as energy 
losses minus energy gain via the radiative processes and q is 
the conductive heat flux. 

The solution for vertical disc/corona structure depends 
on what kind of physical processes are taken into account. 
In the hot corona Comptonization and bremsstrahlung plays 
main role, but if we want to find the structure of the tran- 
sition layer it is obvious that going closer to the cool, dense 
disc also atomic absorption and emission becomes impor- 
tant. Furthermore, the accreting matter heats itself by vis- 
cous forces. Also thermal conduction in vertical direction 
gives substantial contribution to the heating/cooling since 
the transition layer is characterized by strong temperature 
gradient. All physical processes which are included in our 
energy balance equation are described below in section |!| 



3 ENERGY BALANCE - EQUILIBRIUM 
3.1 Radiative processes 

Assuming only radiative energy exchange between the lay- 
ers, the condition of thermal equilibrium is satisfied by: 



£(p,T)=0. 



(8) 



Considering following processes: atomic cooling (spon- 
taneous emission, bremsstrahlung), atomic heating (absorp- 
tion in lines and continuum), Compton heating and cooling 
and viscous heating, the energy equation becomes: 

Ah 

pA(T) - F hard j(T) - Kes F tot 5 (Tic -T) - 



2 



-T = 0, 



(9) 



where A(T) and j(T) are cooling and heating functions re- 
spectively. We use the same A(T) and f(T) functions deter- 
mined by CLOUDY code as in our previous paper (RC96) 
which contains full discussion how to obtain these relations. 
The first function represents the total cooling rate via the 
bremsstrahlung, bound-bound and bound-free emission in 
ergcm 3 s _1 g -2 units, the second one describes the absorp- 
tion in lines and continuum in the cm 2 g -1 units. 

For radiative processes the energy balance equation is 
algebraic and we can expect more than one solution for 
temperature. In fact, we concluded in RC96 that the ver- 
tical structure of the transition zone is not continuous, if 
hydrostatic equilibrium is imposed. For two narrow ranges 
in optical depth the matter is described by three different 
values of temperature and one of them correspond to ther- 
mally unstable solution. Other two stable solutions cannot 
be matched continuously. The transition therefore might be 
discontinuous and located arbitrarily within three-solution 
zone (as suggested by Koo & Kallman 1994) or, as we argue, 
the state of gas is determined by a two-phase equilibrium. 
We concluded there was a narrow slab of matter between 
disc and corona where cool clouds could coexist with hot 
plasma under the constant pressure. 

To describe the nature of the two-phase medium it is 
convenient to use the ionization parameter S defined after 
KMT (1981) as: 



nknTc' 



(10) 



where F ion /c is the pressure of the ionizing radiation, c is 
the velocity of light, n is the number density of particles 
[cm -3 ] and Ub - Boltzman's constant. 

For very large 2 only Compton processes are impor- 
tant and the temperature approaches an asymptotic value 
(i.e. inverse Compton temperature) independent of H. With 
decreasing ionization parameter bremsstrahlung gradually 
becomes more important, leading to the temperature de- 
crease. E£ is the minimum value of ionization parameter 
when high temperature equilibrium is possible. For lower 
S bremsstrahlung cooling overwhelms photoionization heat- 
ing. 

Similar situation is observed for small H. The temper- 
ature first increases with increasing ionization parameter 
(because of increasing of ionization level) until ionization 
parameter reaches critical value H*. At this point photoion- 
ization overwhelms the line cooling and low temperature 
equilibrium becomes unstable. 

For narrow range of ionization parameter H ; * < 5 < H* 
cold, dense matter can exist in pressure equilibrium with 
hot, less dense gas. Both branches are connected by a third 
one with the negative slope which corresponds to thermally 
unstable solutions. 

The ionization state of the irradiated gas in thermal 
equilibrium is described as a relation T(E) which depends 
on: ionizing flux and the shape of cooling and heating curves 
A(T) and 7(T). In our situation hard X-ray radiation is 
predominantly responsible for the ionization of the gas, and 
for the eventual presence of the multi-phase region so we 

accept Fi on — Fhard- 
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3.2 Thermal conduction 

Trying to solve disc/corona transition we should expect large 
temperature gradient so free electrons can efficiently trans- 
port the heat from upper layers to the lower ones. Consid- 
ering such a situation we solve the equation ([?]) where pC 
contain all radiative processes described above. 

Usually the conductive heat transport is treated in two 
limits, depending on the effective mean free path of electrons 
in the medium. The classical thermal conductivity is based 
on the assumption that the mean free path is short in com- 
parison with the temperature scale height T j | VT |. For a 
plasma of cosmic abundance the conductivity is (Draine & 
Giuliani 1984): 



k = 5.6 x 1O" 7 C T 5/2 , 



(11) 



where <f> c is the factor corresponding to reduction in the 
mean free path due to magnetic fields and turbulences (it is 
taken to be <j> c = 1 for equal ion and electron temperature). 
The heat flux in classical case is expressed by the diffusion 
approximation: 



Qclass 



dT 

dz 



(12) 



When the mean free path is comparable to or greater 
than the temperature scale height the heat flux becomes 
'saturated' and nonlocal theory of heat conduction is re- 
quired. In this situation the heat flux depends on the elec- 
tron distribution function and is limited by two constraints. 
First is that divergence of the current must vanish, and sec- 
ond, that electrons streaming through the ions should be 
stable against various plasma instabilities (like ion acoustic 
instability). To estimate of the flux we rewrite after Cowie 
and McKee (1977): 



50s pc = 5(j> s cp, 



(13) 



where c 2 = p/p is the isothermal sound speed, (j>s is the 
dimensionless uncertainty parameter. For Maxwellian dis- 
tribution of electrons cj> s = 1.1, but experimental evidences 
show that usually 4> s ~ 0.3 (Max et al. 1980) and we will 
use this value in our calculations. 

For the purpose of numerical computations, it is conve- 
nient to have the expression which gives smooth transition 
from classical diffusive to the saturated transport. So we 
adopt (after Balbus & McKee 1982) very useful formulae 
defining effective heat flux as: 



k dT 



(14) 



1 + a dz 

where a is the ratio of the classical to the saturated heat 
flux at a given temperature: 



qcla 



Qsat 



(15) 



The energy balance is the differential equation of sec- 
ond order in this case and we expect the temperature to 
be the monotonic function of the optical depth. This means 
the thermal conduction can in principle suppress two-phase 
medium in the transition zone between disc and corona. The 
temperature profile T(t) could be smooth or very sharp. 
Smooth one secures stability of the vertical structure solu- 
tion, but the sharp one should be studied more carefully. 



Perhaps small perturbations of thermodynamic parameters 
can lead to thermal instability and cloud condensation. 



4 LOCAL THERMAL INSTABILITY 
4.1 General considerations 

The thermal stability of the matter interacting with radia- 
tion field was discussed by Field (1965). He assumed a per- 
turbation of thermodynamic parameters in uniform medium 
which in equilibrium state attains L(po,To) = 0. He ne- 
glected any energy transport and dynamic flows, so his iso- 
baric criterion for thermal instability was the following: 



(-) 



<0, 



(16) 



where (d£/dT)p is the isobaric evolution of the loss-heat 
function from the equilibrium. 

Balbus(1986) made more general picture since he took 
into account the dynamical perturbations. In this case the 
loss-heat function of the medium can be a function of den- 
sity, temperature and possibly space and time. Thus the 
isobaric instability criterion becomes: 



/ dln£ \ 
UlnTA 



< 1. 



(17) 



We analyze a local thermal instability taking into ac- 
count the energy transport. Our medium in this case is char- 
acterized by temperature To(z) and density po(z) profile. We 
neglect any dynamic flows and we assume that Tb(z) is at- 
tained by balance between radiative processes and thermal 
conduction. Our unperturbed medium is in stationary state 
denoted by us as U(po,To) and it satisfies the condition: 



p U{po,T ) = p C(p ,T ) + ^ = 0. 



(18) 



When total cooling exceeds total heating (via radiative pro- 
cesses and thermal conduction) then U > 0. We perturb two 
independent thermodynamic variables (density and temper- 
ature) under the condition of constant pressure. We are in- 
terested only in isobaric perturbations because they can lead 
to clouds condensation. The specific entropy of the material 
will change by an amount 5S. This change evolves in time 



d 5S = s d_s = s (g\ = _ § (U 

dt dt \TdtJ \T 



(19) 



If 8S > then its time derivative should be less than zero 
to return perturbed entropy to the background one ( sta- 
bility). For -^SS with positive sign the positive perturba- 
tion will grow up from the background entropy (instability). 
Therefore the instability criterion follows: 



d(U/T) 
dS 



< 0. 



(20) 



Taking into account that in isobaric perturbation 
TdS = CpdT and U = for stationary state, our criterion 
reduces to: 



(dU\ 



(21) 
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4.2 The perturbation growth rate 

In order to find the possibility of the existence of unstable 
regions we perform the local thermal stability analysis. It 
means that we assume our medium to be homogeneous lo- 
cally. This assumption is done, because our problem is local. 
We cannot find full disc structure from the bottom to the 
disc surface for the reason that below our transition layer 
there is another thermally unstable zone as a result of ra- 
diation pressure domination (Pringle, Rees & Pacholczyk 
1973, Shakura & Sunyaev 1976). Global analysis does not 
allow us to separate this two effects. 

At each distance z from the disc surface z we are looking 
for perturbations growing up in time For this purpose we 
consider the perturbation of two independent parameters: 
density and temperature, in the form of a wave in space and 
time: 

T = To + Ti exp(ikz + iut), p = po + pi exp(ikz + iujt), 

where k is a wavenumber and u) is frequency. All quantities 
with index characterize unperturbed state. 

Field (1965) considered such a perturbation of thermo- 
dynamic parameters in homogeneous medium and he de- 
rived the critical wavenumber above which the thermal con- 
duction between perturbed element and surroundings sup- 
presses the thermal (radiative) instability: 



k c — 



k \dTjp. 



(22) 



The inverse of critical wavenumber gives the character- 
istic length scale over which perturbation grows up in a re- 
sult of thermal instability. Let us to denote this length scale 
as A* = l/k c . A A* differs from usually used Field length 
which is defined as: 



Xi 



Tk 



1/2 



(23) 



where Cm = max (cooling, heating) > 0. 

Physically it means that if perturbed element has a ra- 
dius r < A*, then such a cloud will always evaporate due 
to thermal conduction. For perturbation with r > A* cloud 
can condensate depending on radiative conditions in the sur- 
rounding medium (McKee & Begelman 1990). 

For our non-uniform medium with energy transport we 
still assume only isobaric perturbations dP = and we ne- 
glect dynamical flows (so eulerian and lagrangean perturba- 
tions are equal). The matter is characterized by time depen- 
dent energy balance equation: 



3dP 5Pdp T ,,dq 



: 0. 



(24) 



For thermal instability analysis we use the heat flux in clas- 
sical diffusion approximation (equation Jl2f) . 

The first two terms of equation (j24|) are pT(dS/dt) cal- 
culated for perfect gas. 

Considering only linear deviations from the equilibrium 
we derive the relation between ui and k putting perturba- 
tions of density and temperature to equations of energy bal- 
ance and state. The growth rate for each value of z is: 



K fc 2 +po (fr); 
5pok B /m H 



(25) 



Here iuj is always real either negative, or positive. 
We can write temperature deviation as follows: 

ST — Tiexp(Re(iuj)t)exp(ikz) 



(26) 



If the real part of iu> has positive sign then the perturbation 
will tend to increase, for negative sign perturbation will be 
suppressed. Therefore we can find the critical wavenumber 
denoted here by krc which bounds this two cases. From 
condition Re(icu) = 0, krc is given by: 



krc 



po(dU\ 
' k \dTJ P . 



1/2 



(27) 



This relation reduces to the one of Field (1965) in the case 
of homogeneous medium when £o(po,7o) = 0. For perfect 
gas we have: 



fdU\ = (d£\ Co ^fdC 
\dTjp ~ \dTj p To To [dp 



(28) 



The characteristic scale length above which the thermal 
instability develops can be written as Atc = 1/ferc- So we 
expect the perturbed matter with size r > Arc would con- 
dense in clouds with different density and temperature than 
the surrounding medium. We guess that such a perturbation 
can form a two-phase region in the non-uniform matter. 

Local linear stability analysis is valid until exp(Re(iu)i) 
is small enough to neglect second order terms. It gives strong 
constraint on time over which this analysis applies. The time 
scale of linear perturbation cannot be longer than: 



(29) 



Re(iui) 



Computation of further evolution would require global 
stability analysis. This in turn, would be possible only af- 
ter the unperturbated full vertical structure is computed. 
Such complete models are not available at present since ei- 
ther surface layers or deep layers are ignored in the course 
of computing (e.g. Sincell & Krolik 1997, Rozariska et al. 
1998). The results of local stability analysis cannot there- 
fore describe the final state of the transition zone. However 
they may show the existence of the thermal instability zone 
and strongly indicate the possibility of the formation of the 
two-phase medium. 



5 RESULTS 

5.1 Parameters and Boundary conditions 

To solve the structure of the transition layer we integrate 
differential equations presented above as functions of r start- 
ing from the arbitrary (but low enough) values of the optical 



-"gem" 3 . 



depth and density as: to = 2 x 1CP J , po = 10 

In a pure radiative case the temperature is determined 
by solving an algebraic equation (^) for each optical depth. 
In the case with thermal conduction we integrate second 
order differential equation (Q) and the initial temperature 
on the top is computed from equation (^). Second boundary 
condition is set as the thermal flux being zero on the surface 
(dT/dr = 0). 

Soft flux generated via viscous dissipation FdissO, Fhardo 
and X-ray radiation spectrum at the disc surface are free 
parameters in our model. We consider in detail a case of a 
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disc emitting soft flux Fduao = 5 x 10 ergcm s which 
corresponds to the black body temperature 3.16 x 10 4 K. 
Such a flux is expected from a disc around massive (mass of 
10 8 Afo) black hole at WRschw if the accretion rate within a 
disc itself is about 2.8 x 10 24 g/s (i.e. about 0.01 in Eddington 
units) . 

We assumed X-ray radiation spectrum with energy 
cutoff at 200keV, as usually observed in Seyfert galaxies 
(Zdziarski et al. 1995) and with spectral indexes: oe = 
-0.987 for hv < 200keV and a E = -2.1 for hu > 200 keV. 
Results are presented for different values of hard radiation 
flux: F hard o = 5 x 10 14 and 5 x 10 15 ergcm" 2 s _1 . 

In the case without thermal conduction the thickness 
of the accretion disc was adopted as a minimum possible 
value H drain for which the disc is optically thick (RC96). 
For Hd < H 'dmin disc becomes transparent, whereas for 
Hd > Hdmin disc is no more geometrically thin and hot 
atmosphere becomes very narrow in comparison with cold 
matter. Therefore we adopt the disc thickness Hdmin in our 
approach. 

Other quantities adopted in our calculations are typ- 
ical for standard accretion disc model. We take the vis- 
cosity parameter equal a = 0.1, the Keplerian angular 
velocity is SI = 10~ 



5 s 1 (this value corresponds to Kep- 



hole at the mass of 10 Mg), opacity coefficient for elec- 
tron scattering is K es = 0.34 cm 2 g _1 and Kramer's co- 
efficient Ko = 3.8 x 10 23 cm 2 g _1 . Total opacity coefficient 
used in relation ^| is computed from the standard relation 

Ktot — K es + KopT ' . 

The results are presented as the functions of optical 
depth t or of the distance from the disc surface z. Both 
parameters are connected by relation: 



dz 



(It 



Ktot(T)p' 



(30) 



5.2 The vertical structure of the transition layer 

The computations show that in our model heat flux is carried 
according to classical description, because a <C 1 (see eq. 

§)• 

Cool, optically thick disc is almost isothermal so, the 
thermal conduction between the layers is negligible. Simi- 
lar situation is observed in hot corona, but some heat flux 
can appear there until corona has really hight temperature 
(classical conductivity is proportional to the T 5 ^ 3 ). Our nu- 
merical results show that heat flux generated in corona is 
about six orders of magnitude smaller than the one in tran- 
sition layer. 

Fig.l. presents the temperature profiles T(r) (I) and 
stability curves T(E) (II) of the disc illuminated by: a) 
Fhardo = 5 x 10 14 ergcm" 2 s _1 and b) F hard0 = 5 x 10 15 
ergcm~ 2 s -1 . Computations for the case when thermal con- 
duction is taken into account (solid line) are compared to the 
ones without thermal conduction (points). Near the surface, 
the disc atmosphere for both cases is heated up to the same 
temperature for given Fhardo- The value of this temperature 
is close to the Inverse Compton temperature which is de- 
termined by the spectral shape of incident X-ray radiation. 
In our case Tic — 8.19 x 10 7 K on the top of the discussed 
layer. The zone with hot plasma is thicker for larger F nar do- 
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Figure 1. Temperature versus optical depth r (left hand side I) 
and ionization parameter H (right hand side II) for two differ- 
ent values of incident hard radiation flux: a) F na rdO = 5 X 10 14 
ergcm - 2 s — b) Fhardo = 5 X 10 15 ergcm - 2 s _1 . Solid wide line 
represents the case with conductivity, thin dotted line shows the 
pure radiative case. 



When the density increases with the optical depth the 
importance of the bremsstrahlung cooling increases and the 
temperature starts to fall down. For the case with conductiv- 
ity and for lower Fhardo this point appears at lower optical 
depth than for pure radiative case. 

It is shown in previous paper (Figure 5 of RC96) that 
our method is applicable for r > 1 and we do not have to 
consider the radiative transport equation more carefully. 

Going further inside the disc temperature reaches the 
intermediate stable region at T ~ 8 x 10 5 K. Since photoion- 
ization is the predominant process in the thin transition 
layer between hot T ~ 10 8 K and cold T ~ 10 4 K matter, 
we can expect that heavy elements abundance influences 
the shape of instability curve. For some configuration of ele- 
mental abundance even two intermediate stable regions can 
exist (Hess et al.1997). In our model one intermediate equi- 
librium region forms and it is wider for larger values of X-ray 
radiation flux. 

The most interesting feature of the presented profiles is 
that for the pure radiative case there are two narrow ranges 
of r when matter can reach three different values of temper- 
ature and one of them is unstable dT/<E < 0. For the case 
with conductivity the profiles are very sharp but continuous 
opening a possibility that the thermal conduction suppresses 
thermal instability in the transition layer between disc and 
coronae. 

Finally, the temperature drops to the value of 2.8 x 10 4 
characteristic for the cool disc. Bigger flux of X-ray radia- 
tion illuminating the surface leads to higher optical depth 
reached by cold matter. Our computations for the case with 
conductivity are stopped at the point where the thermal 
conduction flux vanishes at low temperature equilibrium. 

Figures named by II (right hand side) show ionization 
curves and we can easily see that thermal conduction does 
not change the ionization state of matter in regions of ther- 
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Figure 2. Vertical distribution of I) pressure, II) density and 
III) optical depth at two different values of incident hard radia- 
tion flux: left hand side a) F hard Q = 5 X 10 14 ergcm - 2 s — 1 , right 
hand side b) ffc OP( jo = 5 X 10 15 ergcm _2 s — 1 . Solid wide line rep- 
resents the case with conductivity, thin dotted line shows the pure 
radiative case. 



mal stability. But two discontinuous parts of these lines 
which are characteristic for the pure radiative case are re- 
placed by ones with positive slope dT/dS > 0. So the solu- 
tion of the disc structure for stationary state is very steep, 
but monotonic. 

Fig. 2. presents variations of other disc parameters like: 
pressure (I), density (II), and distance from the surface z 
(III) with optical depth. We again compare pure radiative 
case (points) and the case with conductivity (line) for two 
values of F^ardo a) and b) which are described above. 

Near the surface the profiles are identical for both cases. 
At some optical depth (different for each parameter) lines 
start differ and when thermal conduction is taken into ac- 
count gas attains higher pressure and density at lower dis- 
tance from the surface than without the thermal conduction. 
This difference is bigger for lower value of the initial hard 
radiation flux a). 



5.3 Thermal instability criteria for the case with 
conductivity 

Now we would like to focus on the case with thermal con- 
duction, because it seems to be more realistic. It is obvious 
that in medium with strong temperature gradient conduc- 
tive heat transport plays a big role and influences the struc- 
ture of region. In our model there are two very narrow (first 
dr = 1 x 10 -4 and second dr = 1.6 x 10~ 6 ) zones with sharp 
temperature and density profiles in the transition layer be- 



tween an accretion disc and a corona. Thermal heat flux 
rises to 3 x 10 8 in these places. 

To discuss the stability of such a solution we use local 
theory developed in Section 4. As was mentioned in above 
this theory is valid for short period of time comparable with 
1 /Re(iuj) . Otherwise we should treat our problem more glob- 
ally which is more complicated and is not presented in this 
paper. 

Another assumption of applicability of local theory is 
that perturbation of considered parameter x is smaller then 
characteristic scale length xj \ Vx | on which this parame- 
ter changes. Therefore we are looking for any perturbations 
bigger than Xtc, but smaller than T / | VT |, because they 
can lead matter to condensate into clouds. Treating prob- 
lem locally we determine the dependence of Xtc and the 
temperature scale length (named here Ht) on the optical 
depth. 

For a wide range of r the solution is stable, i. e. k TC < 
(equation (p7|)). Only for those two narrow zones with ex- 
tremally strong temperature gradient, wavenumber becomes 
positive and we can find size of perturbations that lead to 
thermal instability. Fig. 3. presents variation of Xtc and Ht 
versus r for first (I) and second (II) zone, and for two dif- 
ferent values of hard radiation flux Fhardo — 5 x 10 14 a) 
and F hard0 = 5 x 10 15 b). Because possible unstable regions 
are very thin, we rescale drawings by subtracting the optical 
depth of the beginning of zones to from r. Values of to are 
written in the lower left corner for each case. Note that the 
optical thickness dr of zones does not depends on X-ray ra- 
diation flux. Only optical depth r and distance of the zones 
from the disc surface z is bigger for higher F^ardo- 

Characteristic temperature scale length is always one 
order of magnitude bigger than Xtc for the first zone in 
case al) and more than half of order of magnitude in case 
bl). So the perturbations with size A where Xtc < X < Ht 
can lead to cloud condensation. The half order of magnitude 
difference between Xtc and Ht is observed for the second 
zone, but not from the beginning where there is no place 
even for minimal perturbations (all, bll). The higher hard 
radiation flux the smaller perturbations are required to pro- 
duce thermal instability in the transition layer. 

We expect our clouds to be under influence of gravity 
and particularly radiation pressure. The life time of clouds 
is very short and probably they don't have enough time to 
grow up to the much bigger size than the minimal size of 
perturbations. So the most probable size of clouds is in order 
of Xtc and usually is taken into account in calculations of 
cloudy regions (Kurpiewski et al.1997, Krolik 1998). As we 
can see in unstable zones there is not space for a big clouds. 
But this is still enough for several clouds in size comparable 
with Xtc and because they are most important the results 
give good evidence for instability existence in the transition 
layer. 



6 DISCUSSION 

In this paper we described the influence of thermal con- 
duction on the vertical structure of the transition layer be- 
tween a hot corona and a cool accretion disc. We found that 
heat transport via conductivity does not suppresses thermal 
instability which occurs in the case when energy exchange 
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Figure 3. Characteristic length of perturbation Xtc versus opti- 
cal depth (solid line) compared with temperature scale length Ht 
(dotted line) for first (I) and second (II) unstable zone for two dif- 
ferent values of incident hard radiation flux: a) FhardO = 5 X 10 14 
ergcm^ 2 s _1 , b) F^ardO = 5 X 10 15 ergcm _2 s -1 . to is the optical 
depth of the beginning of zones. 

through the disc/corona boundary is only due to radiation 
processes. 

The temperature profile for the case with conductivity 
is continuous but very sharp particularly in the place where 
thermal instability was found for pure radiative case. Local 
stability analysis of stationary state with thermal conduc- 
tion leads to the conclusion that in zones with extremely 
sharp temperature gradient, the spontaneous perturbations 
with size A, where Xtc < A < Ht can cause cloud conden- 
sation. The estimated size of possible clouds is ~ A and is 
much smaller than the size of thermally unstable zone in 
pure radiative case which occurs ~ 10 10 cm for both cases of 

FhardO ■ 

Therefore we expect the structure from the corona sur- 
face to the disc midplane as follows: first the outermost zone, 
hot corona cooled by inverse Compton emission (Begelman, 
Mckee & Shields 1983, Begelman & Mckee 1983, Ostricker, 
Mckee & Klein 1991), second, the transition layer with possi- 
bility of the development of the two-phase medium discussed 
in this paper, and third, main body of the disc (Shakura & 
Sunnayev 1973). We argue that upper parts of the transition 
consists of cool clouds at column density Nh > 10 22 embed- 
ded in hotter medium, but closer to the disc the situation 
changes continuously in the opposite one when hot blobs at 
Nh ~ 10 19 are deep in cooler matter. 

Using local stability analysis we are not able to find the 
evolution of two-phase zones. We can only find the possibil- 
ity of clouds existence. We cannot to integrate full vertical 
disc structure, because ist is very difficult to describe all 
physical processes, particularly in the disc interior (three 
body recombination). Also on the distance from the Black 
Hole equal Wrsch disc is radiation pressure dominated and 
another thermal instability plays role. Global analysis does 
not allow us to separate these two effects . (In the disc in- 
terior convection may be important as well.) In addition we 



found that this is not Sturm-Liouville problem, because we 
can find only finite number of eigenvalues (this is against 
main assumption of Sturm-Liouville problem) . In the future 
we should integrate equations numerically without any as- 
sumptions about character of perturbations and about dif- 
ferential operator. 

These results should be treated as preliminary and fur- 
ther research is necessary first to improve the physical con- 
tent of the model and and second, to explore the range of 
parameters. In the future work the transfer of hard radiation 
through the disc matter should be consider more carefully 
especially in the aim of better determination of absorption 
and reflection of X-rays depending on optical depth. 

The position of possible unstable regions depends on 
the shape of atomic cooling and heating rates. We tested the 
accuracy of the adopted functions in our previous paper. Re- 
cent studies on photoionization codes show that abundance 
of iron and oxygen ions (with smaller contribution of magne- 
sium, silicon and sulphur) plays crucial role in the nature of 
instability. Varying of this abundance strongly influences the 
heating and cooling behavior and the suppression of thermal 
instability might be achieved through small changes (Hess 
et al.1997). 

Our results are important in the aim of determine 
strictly the position of cloudy boundary layer. According to 
the present results, possible unstable regions in stationary 
state occur for 3 = 0.28 (first zone) and 3 = —0.19 (second 
zone). Improved models should be used in the future to de- 
termine the dependence of these values on accretion disc pa- 
rameters (including radius), X-ray radiation spectrum and 
chemical composition of the gas. The global stability analy- 
sis should be done by solving second order differention equa- 
tion and finding global modes of thermal instability. 

We do not consider dissipative corona above transition 
layer. The heat transport via conductivity from such hot 
~ 10 9 K region can change our picture. We expect to do it 
in the future. 

Unstable zones might be very dynamical and it is dif- 
ficult to study the time evolution of such a region. Clouds 
can evaporate or condense or to be accelerated by radiative 
pressure. The first attempt to compute the spectrum from 
cloudy transition layer was done by Kurpiewski et al. (1997). 
They found that broad emission lines observed in quasars 
spectra can come from clouds which form continuously at 
the basis of the corona because of thermal instability. 
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